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Abstract 



Numerical simulations of the two-dimensional t- J model in the limit J/t<l 
are performed for rather large systems (up to N = 12 x 12) using a world- 
line loop-algorithm. It is shown that in the one-hole case with J = 0, where 
no minus signs appear, very low temperatures (fit ~ 3000) are necessary in 
order to reach Nagaoka's state. J/t^0.05 leads to the formation of partially 
polarized systems, whereas J/t>0.05 corresponds to minimal spin. The two- 
hole case shows enhanced total spin up to the lowest attainable temperatures 
{fit = 150). 

PACS numbers: 75.10.Jm 
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Since the discovery of high-T c superconductors (HTS), a great deal of interest was focused 
on strongly correlated systems like the Hubbard and t-J models. Especially the latter 
system attracts a lot of attention since it can be viewed as a model on its own right for 
HTS [0 and not only as the strong coupling limit of the Hubbard model. A key method 
for their understanding is provided by numerical simulations. For the Hubbard model, 
quantum Monte Carlo (QMC) simulations [[| can be performed for not too large values 
of the interaction U 12 in units of the hopping matrix element t) and also not too 
large inverse temperatures (3 (<10/t), except for the half-filled case, where the so-called 
minus-sign problem is absent, limiting thus our understanding of the doped phase at low 
temperatures in the strong coupling limit. On the other hand, exact diagonalizations of the 
t-J model || provide information of precisely that limit but confined to rather small sizes N 
(^S 30). QMC simulations for the t-J model were successfully performed in one-dimension [|J, 
whereas previous simulations in higher dimensions suffered under metastability problems || 
or were restricted to the computation of energies in the ground state ||. Further numerical 
studies were carried out using high-temperature expansions (HTE) |7|,§] and density-matrix 
renormalizat ion-group (DMRG) |9|,|I0]|. 

In this Letter QMC simulations of the t-J model are presented for the first time for large 
systems (N < 12 x 12) and a large temperature range. Large system sizes are necessary 



especially in the case J = 0, as shown by Barbieri et al. [11], where a change in the stability 
of the fully polarized ferromagnetic state is observed for two holes at around N = 12 x 12. 
The simulations are based on a recently developed world-line loop-algorithm | 12| . We use 



a representation of the t-J model where it is obvious that it has no minus sign problem 
in the case of one hole and J = ||i3|,|T4||. As shown below, this enables us to reach very 



low temperatures ((3t ~ 3000), that are necessary to converge to Nagaoka's ferromagnetic 
ground state. Such a slow convergence suggests that a high density of low lying excitations 
is present. We interpret them in the frame of a spin- wave theory, such that an effective 
ferromagnetic exchange interaction J e ff can be assigned for each lattice size. Finite-size 
scaling shows, that J e ff — > in the thermodynamic limit as is expected. Hence, Nagaoka's 
ferromagnetic state has zero spin-stiffness in that limit. A frequently addressed question is 
the stability of Nagaoka's state for J ^ or more than one hole [|T5|-|TT||. Our simulations 
show that for N = 10 x 10, J > 5 x 10 -4 suffices to bring the system to a ground-state without 
maximal spin. Due to the minus-sign problem, simulations for more than one-hole could not 
be performed at temperatures as low as in the one-hole case (j3t < 150) and less conclusive 
results could be obtained. However, it should be remarked that such temperatures are by far 
lower than those reached previously by either QMC (fit <2) || or HTE <5) ||. Finally, 
on increasing J, we show that the system goes to minimal spin at J /£> 0.05. At smaller J 
the system is partially polarized. 

We first describe shortly the representation of the t-J model used in the simulations. 
Exact operator identities lead to a decomposition of the standard creation (c| s ) and anni- 
hilation (cj jS ) operators for fermions with spin s =f, I |]T3| , |l~4 



c !,t = li,+fi ~ 7i -fi » 4,1 = <Ti-(fi + fi) > W 

where 7 ij ± = (1 ± cr i)2 )/2, a i ± = (a iiX ± zcrj j2/ )/2, the spinless fermion operators fulfill 
the canonical anticommutation relations {fj,fj} = and a^ a , a = x,y, or z are the 
Pauli matrices. The constraint to avoid doubly occupied states reduces in this case to 
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J2i1i,-fifi — which is a holonomic constraint in contrast to the one normally used with 
the standard representation. Moreover, this constraint commutes with the Hamiltonian (0), 
such that once the simulation is prepared in that subspace, it remains there in the course 
of the evolution in imaginary time. The t-J Hamiltonian in this new representation has the 
following form 

Ht-j = t £ (/'„/;./• + ^{Pa - 1)) , (2) 

<i,j> V Z J 

where — (1 + 0{ • <Jj)/2, Ay = (1 — rij — rij), and rii = The constant added in (H) 
ensures that H t _j reduces to the standard t-J model. In the case of one hole with J = 0, 
there is only one fermion present and the simulation has no minus-sign problem due to the 
exchange of fermions. For J > a minus sign can occur because the exchange of two spins 
in space time generates them. 

Next we describe the world-line loop-algorithm. The world-lines are as usually defined 



by checkerboarding space-time [0. However, a direct application of this algorithm in two 
spatial dimensions is not effective in dealing with identical particles, since each permutation 
of any two particles represents a different sector that cannot be reached from another one. 
Especially in the case of fermions, the minus signs would lead to an extremely innacurate 



estimate of any observable. This difficulty can be overcome with a loop algorithm [12 
which does not preserve the linking topology of the world-lines. Since we are dealing with 
three states per site, two world lines are introduced (one for the fermions and one for the 
spin-ups) and their updating is performed separately, the sequence being chosen at random. 
Technical details will be published elsewhere (TP] . 

In the following we describe the results for the limit J — > and mainly with one hole. 
Figure [I] shows the structure form factor 




S (t) = E \™P (»*« ■ 9) E SjSfa (3) 



at q — 0, i.e. < S^ ot 2 > for N = 8 x 8 and N = 10 x 10 sites, as a function of fit. As 
can be clearly seen, very low temperatures are needed in order to reach the fully saturated 
ferromagnet as predicted by Nagaoka's theorem. This suggests that a high density of low 
energy excitations is present as would be expected in the case of a Heisenberg ferromagnet. 
As a simple test of this hypothesis, we compare the temperature dependence of S(q = 0) from 
the QMC simulation with a spin- wave result, where the ferromagnetic exchange coupling 
J e ff of the hypothetical Heisenberg ferromagnet is chosen such as to fit the results of the 
simulation in Fig. |l[ Figure |2| shows that J e ff ex 1/N, such that in the thermodynamic limit 
the spin-stiffness vanishes. 

Figure [3] shows that, at the attainable temperatures, the ferromagnetic correlations are 
enhanced by adding a second hole to the system (N = 12 x 12), suggesting a picture of a 
ferromagnetic polaron surrounding each hole, with a size limited by thermal fluctuations. 



Our result is consistent with a stability analysis by Barbieri et al. ['II , where the Nagaoka 



state is stable for two holes and 12 x 12. However, such a comparison should be taken 



with care since the results of Ref. [11] are strictly valid in the limit \nN, where 



is the number of holes. Also DMRG studies |l0j support a Nagaoka state for such a low 
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doping (5 ~ 1.4%), although the geometry used in that case is quite different from ours. 
Unfortunately, the minus-sign problem precludes simulations at lower temperatures, such 
that the ground state cannot be reached. In fact, by comparing Fig. |3| with Fig. [I], it is clear 
that in order to extract conclusions about the total spin in the ground-state, temperatures 
much lower than (3t = 150 are necessary, such that previous assertions on the character of 
the ground state from results at [3t = 2 to 5 j5]|| are doubtful. 

Exact diagonalization studies have shown, that the Nagaoka state breaks down for J > 



Ji ~ 0.075t in a 4 x 4 system ||20|| . Further they show, that the system goes to minimal spin 
at J 2 ^ 0.088. Between J\ and J 2 the system is partially polarized. In the following we show, 
that for larger systems (N = 10 x 10) J\ and J 2 are of different order of magnitude. Whereas 
Fig. [| shows that the system is partially polarized for J = 0.0005t, Fig. | demonstrates 
that for J>0.05t the system goes to a state with minimal total spin. At the same time, 
antiferromagnetic correlations begin to dominate at low temperatures. Although the minus 
sign problem does not allow a very accurate determination, our data show, that the system 
is partially polarized for 0.0005 < J < 0.05. The last boundary is somewhat higher than the 
lower bound proposed for the Hubbard-model to reach an antiferromagnetic state in the 
one-hole case (Uf F > (At/n)N\nN) [0. 

Summarizing, we have performed QMC simulations of fermions in the limit of infinitely 
strong correlations. The world-line loop-algorithm allows simulations of rather large systems 
(N < 12 x 12) and very low temperatures (up to fit ~ 3000). Such temperatures are necessary 
in order to reach the fully polarized state in the case of one-hole and J = 0. J/t< 0.05 leads 
to the formation of partially polarized systems, whereas J/t^0.05 corresponds to minimal 
spin. The two-hole case shows enhanced total spin up to the lowest attainable temperatures 
(f3t = 150). 

We are very gratefull to Antimo Angelucci for instructive discussions on the initial stages 
of the project. The calculations were performed at LRZ Miinchen and HLR-Stuttgart. We 
thank the above institutions for their support. 
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FIGURES 



FIG. 1. S(q = 0)/N for N = 8 x 8 and 10 x 10, with one hole at J = 0. The dotted line 
represents the fully saturated ground state. The full lines correspond to S(q = 0) of a Heisenberg 
ferromagnet in spin wave theory with J = 0.0065t and J = 0.0030i, respectively. 



FIG. 2. J e ff of the Nagaoka state for different lattice sizes N. 



FIG. 3. Comparison of one hole and two holes in a 12 x 12 lattice. One sees, that ferromag- 
netism is enhanced by adding another hole. 



FIG. 4. S(q = 0) for a 10 x 10 lattice with one hole at t > J > 0. 



FIG. 5. S(q) for q = and q = (it, ir) for a 10 x 10 lattice with one hole at J = 0.05£ 
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